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We present a lattice computation of the first moment of the kaon's leading-twist distribution am- 
plitude. We use ensembles with 2+1 dynamical flavours of domain wall fermions and the Iwasaki 
gauge action from the RBC and UKQCD joint dataset. We observe the expected chiral behaviour 
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1. Introduction 



We present a lattice calculation of the first moment of the leading-twist distribution amplitude of the 
kaon, 0jf(«,jU) Among the many phenomenological applications which require knowledge of 
distribution amplitudes are electromagnetic form-factors at large momentum transfer and related 
processes |||, ||, |5|, ^, [7l |h, and, following the development of the factorization framework, 
exclusive charmless two-body B-decays into two light mesons U [T(], [II], (T^J, [l3|, p"4, 15]. 
The distribution amplitude parametrizes the overlap of a kaon with longitudinal momentum p with 
the lowest Fock state consisting of a quark and an anti-quark carrying the momentum fractions 
up and Up = (1 — u)p, respectively (u + u = 1). It is defined by the non-local (light-cone) matrix 
element 



(0\q(z)Y P Y5^(z,-z)s(-z)\K(p))\ z2 



z 2 =0 



f K (ip P ) [ due^r^Kiu^), (1.1) 
Jo 



where /I is a renormalization scale and £P(z,—z) = exp {— ig J^ z dw >1 A ll (w)} . The distribu- 
tion amplitude is normalized by J du§K{u,[i) = 1 and can be expanded in terms of Gegenbauer 

3 /2 

polynomials C n (2u — 1), 



§ K (u,pL) =6uu ( 1 + ^ a n {n)C?/ 2 (2u- 1 



(1.2) 



n>l 



The lowest Gegenbauer moment af is proportional to the average difference of the longitudinal 
quark and anti-quark momenta of the lowest Fock state, 

«f(M) = | C du{2u-\)^ K (u,pL) = \(2u-\) = 5 -{^){pL). (1.3) 
3 Jo 3 3 



While the first moment of the distribution amplitude vanishes in the case of the pion, it is non-zero 
for the Kaon because of SU(3)-breaking effects. ( £ ) is obtained from the matrix element of a local 
operator, 



(0\q(0)Y P Y5D^(0)\K(p)) = (^)f K p pPfl = - aff K p p Pfi . 



(1.4) 



d n + igA^ and D M 



igA 



where we use = —D^, 
The first moment of the kaon's distribution amplitude has in the past been determined mainly from 
QCD sum rules, and recent results include: af (lGeV) = 0.05(2) 0.10(12) @], 0.050(25) 
[18] and 0.06(3)[|T9|]. Very recently an independent lattice study of this quantity was published [ |2"(]| ] 
which quotes af (2GeV) = 0.0453 ± 0.0009 ± 0.0029 as the final result. 

Here we use the N f = 2+l gauge field ensembles from the RBC and UKQCD dataset [[H], [22], |23|] 



(domain wall fermions [24, 25] and Iwasaki gauge action [26, 27]) with three values of the light- 
quark mass with m sea = m va i ence in each case. The hadronic spectrum and other properties of these 



configurations have been presented at this conference [21, 22, 23] 
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2. ( E, ) bare from Lattice Correlation Functions 

In constructing the lattice operators which are relevant for the determination of ( £ ) , we use the 
following symmetric left- and right-acting covariant derivatives: 

DnY(x) = — W(x,x + fi)\ir(x+fi)-U(x,x-fl)\ir(x-fi)} , (2.1) 
2a 

y{x)D u = —{xjr(x+fX)U(x+fl,x)-Xff(x-fl)U(x-fX,x)} , (2.2) 
2a 

where the IPs are the gauge links and fi is a vector of length a in the direction /J. (a denotes the 
lattice spacing). 

To illustrate the method, consider the local lattice operators O p ^(x) = q(x)YpYs^u s ( x )^ A p (x) = 
q(x) YpY=,s(x) and P{x) = q(x) y$s(x) from which we define the two-point correlation functions 

Cpv(t,P) =I^(0|O pM (f,x)P+(0)|0) and C AvP (t,p) =^e i ^(0\A v (t,x)P^0)\0) . (2 . 3) 

X X 

Here q and s represent the light and strange quark fields, respectively. At large Euclidean times t 
and T — t (T is the length of the lattice in the time direction), we expect 



, f .• , - C {PH}(^P) _^ ; PpPn i e \baie n 

R{p ^ p) = ^m • (2 - 4) 



The superscript bare denotes the fact that the operators are the bare ones in the lattice theory with 
ultraviolet cut-off a~ 1 in the Domain Wall Formalism and the braces in the subscripts {ppi} indicate 
that the indices are symmetrized. In order to avoid mixing of 0^ v under renormalization [^] we 
only consider the cases p = V =4, \i = k (k = 1,2,3) with p% = ±2n/L while \p\ = 2n/L. 



3. Perturbative Renormalization of the Lattice Operators 



The perturbative matching from the lattice to the MS scheme is performed by comparing one- 
loop calculations of the two-point Green function with an insertion of the operator 0{p^} in both 



schemes. Defining 0^3, Li (aO = Zq, , Of^ ^ (a), the renormalization factor is given by 



Zo 



1 - wi)Z u 



l + ^(-^HLL 2 a 2 )+Zr-L l+ V MS -V 



(3.1) 



In this expression, (1 — Wq)Z„ is a characteristic normalization factor for the physical quark fields 
in the domain wall formalism. It is a common factor in the numerator and denominator of the ratio 
^{pji};v as we the contributions from the wave function renormalization. Z w represents an additive 
renormalization of the large Dirac mass or domain wall height M = 1 — wq which can be rewritten 
in multiplicative form at one-loop as Z w = 1 + j^z w with z w = 

The terms Z^ s and £i come from quark wave function renormalization. The terms V MS and V 
come from the one-loop corrections to the amputated two-point function. Using naive dimen- 
sional regularisation in Feynman gauge with a gluon mass infrared (IR) regulator, £^ s = I and 
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MS 



25 
18- 



V 

The contribution Ei has been evaluated for domain wall fermions with the Iwasaki gluon action 
in Feynman gauge in p9||. We have calculated the lattice vertex term V for the same action and 
gauge regulator to complete the evaluation of Zo {pfl} . The perturbative calculation is explained 



in [30, 



j, [3_lJ and the form of the Iwasaki gluon propagator can be found in \p2\]. 
For the Iwasaki gluon action and for the value of M = 1.8 used here the physical quark normaliza- 
tion z.w has been found to be very large in [ |30| , f^] and we therefore use mean field improvement as 
described in [[29|]. 

The first step is to define a mean- field value for the domain wall height, M MF = M - 4(1 - P 1 / 4 ) 
where P = 0.58813(4) is the average plaquette in our simulations, leading to M MF = 1.3029. 



The physical quark normalization factor becomes [l — (w^ F ) 2 ] Z^ F , with Z^ F = 1 



g C F „MF 



,MF 



2w ^- (L w + 32k 2 T, 



1671- 



and 



MFj 



is a mean-field tadpole 



5.2509, where T MF = 0.0525664 

3.9731 and V = -4.1907 in equation (|3T|) are 



factor and L w is evaluated at M MF . Likewise, £i 
evaluated atM MF and the mean-field improved renormalization factor for our simulations becomes: 



1 



0.9082 



1 - 1—4 5.2509 
16k 2 



1 + 



g 2 C F 
16k 2 



-InQuV) -0.6713 



(3.2) 



We make two choices for the mean-field improved MS coupling. The first uses the measured 
plaquette value, P, according to [|29|] 



1 



P 22 

7 2+dg+Cp + T6K~ 2 



(3.3) 

6/g 2 = 2.13 in our sim- 
= 1, we find aj^§-(plaq) = 

0.1752 and aj^g-(ctm) = 0.3385. With these two choices of coupling, our value for the renormal- 



where d g = 0.1053 and c p 



0. 1401 for the Iwasaki gauge action and j3 - 
ulations. The second choice is the usual continuum MS coupling. At \ia 



*MS 

ization factor becomes 



Z A 



1 .2346 plaquette coupling 
1.3384 continuum MS. 



(3.4) 



We include the spread of results in eq.(3.4) as the estimate of our current systematic uncertainty in 

— 1.28 ± 0.05 for the final result. 



the renormalization factor and thus we will eventually use — ^ 



4. Numerical Simulation and Results 

The lattice volume is (L/a) 3 xT ja x L s = 16 3 x 32 x 16. The choice of bare parameters is j8 = 2.13 
for the gauge coupling, am s = 0.04 for the strange quark mass (which has been tuned to correspond 
to the physical value) and am q = 0.03, 0.02, 0.01 for the light-quark masses. With these simulation 
parameters the lattice spacing is cT l = 1.60(3) GeV [|22|, ^3|]. Owing to the remnant chiral symme- 
try breaking the quark mass has to be corrected additively by the residual mass in the chiral limit, 



am res = 0.00308(3) [22, 23]. 
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Figure 1: Left: Jack-knife results for ( | } bare as a function of the time. The ranges over which we fit and the 
corresponding results are indicated by the black lines. Right: Linear chiral extrapolation for ( t, } bare . 

4.1 Bare correlation functions 

For each value of the light-quark mass we computed the correlation functions on 300 gauge config- 
urations separated by 10 trajectories in the Monte Carlo history. On each configuration we average 
the results obtained from 4 positions of the source for the lightest quark mass (am q = 0.01) and 2 
positions of the source for the remaining two masses (am q = 0.02 and 0.03). In order to improve 
the overlap with the ground state at the source where we insert the density P*, we employed gauge 



invariant Jacobi smearing [ |33| ] (radius 4 and 40 iterations) with APE-smeared links in the covariant 

0.4164(10), 



Laplacian operator (4 steps and smearing factor 2) [J34j, p5| ]. 

The kaon masses corresponding to the simulated bare light-quark masses are ara^ 03 
am° K 02 = 0.3854(10), and am° K 01 = 0.3549(14). 

The left plot in figure [T] shows our results for ( t, ) bare as a function of t obtained from the ratio 
R{Ak}\A^iPk = ±27r/L) for the three values of the mass of the light quark. We averaged the results 
over equivalent choices for the momenta and folded the data in the time-direction. There are clear 
plateaus, demonstrating that the SU (3) -breaking effects are measurable and ( E, ) can be determined. 



4.2 Chiral extrapolation 

Plotting our results for ( E, ) bare as a function of the light-quark mass in the right plot in fig. [I] and 
taking into account the remnant chiral symmetry breaking by defining the chiral limit at the point 
am q + am Tes = our data confirms the linear behaviour predicted by chiral perturbation theory p6| , 



37]. Moreover the line passes through ) bare = at a value of the light-quark mass (denoted by 
the open square) which is consistent with the mass of the strange quark, as expected for the SU (3) 
symmetric case (am uc i = am s = 0.04). From the linear fit we obtain ) bare = 0.0262(23) in the 
chiral limit. 

5. Systematic Uncertainties and our Final Result 

Combining ) bare with the result for the perturbative renormalization factor we obtain our final 
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result 

(£) m (H = 1. 6 GeV) =0.034 ±0.003. (5.1) 

In order to compare our result with previous calculations we evolve it to the renormalization scales 
1 GeV and 2 GeV using the three-loop anomalous dimension [ 38 ] . We obtain ( t, ) MS (fi = 2 GeV ) = 
0.032 ± 0.003 and ( % ) m (n = 1 GeV) = 0.040 ± 0.004. 

The error in the renormalization factor due to the uncertainty in the lattice spacing is negligible. 
For example if we conservatively allow the lattice spacing to vary between 1.58 GeV and 1.62 GeV, 
the contribution to the relative error on ( £ ) MS is less than 0.2%. 

Among the uncertainties which we are not at this stage in a position to check numerically are the 
continuum extrapolation, finite-volume effects and the fact that the strange quark mass (m s a = 0.04) 
is only approximately tuned to its physical value. The lattice artefacts are formally of 0{a 2 A^ CD ) ~ 
2.5% and we are planning to check this with a simulation at a smaller lattice spacing. We would 
expect the finite volume effects to be small and are currently checking this with a simulation on 



a 24 3 x 64 lattice. The strange quark mass appears to be well tuned [22, |||] so again we expect 
the contribution to the error from this uncertainty to be very small. Thus we expect the errors 
from these three sources to be sufficiently small not to change the errors quoted for our final result. 
We are also carrying out a systematic programme of non-perturbative renormalization which will 
enable us to reduce the uncertainty in the renormalization constants. 



6. Summary and Conclusions 



We have demonstrated that the SU (3) -breaking effects which lead to a non-zero value for the first 
moment of the kaon's distribution amplitude are sufficiently large to be calculable in lattice simula- 
tions and satisfy the expected chiral behaviour. As our best result we quote ( E, ) MS {pL = 1 .6GeV) = 
0.034 ±0.003. 
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